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itly constructed and their relation to separation of variables in H-J 
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1 Introduction 



Few years ago Braden [I] classified all translationally invariant systems of 
identical particles admitting an additional (apart from energy and total mo- 
mentum) integral of motion, which is polynomial in momenta of degree three 
(actually, an additional assumption in Braden's paper is that the third-order 
term is q - independent). It appeared that the assumption of SV symmetry 
put the severe restrictions on the admissible hamiltonians. 
In the present paper we classify the dynamical systems of identical particles 
possesing the quadratic integrals of motion. We find the general form of po- 
tentials and we show that the additional integrals arise from the separation 
of radial variable in Jacobi coordinates. We explain also the group - theoretic 
origin of the integrals of motion. 

2 Quadratic integrals of motion 

We want to classify all classical dynamical systems of iV degrees of freedom 
obeying the following assumptions: 

(i) the hamiltonian takes the natural form 

1 N 

ff=«£F? + V(fc. ■••,3*0 (!) 
/ i=i 

(ii) H is translationally invariant 

{H, P} = 0, P = f> (2) 

i=i 

(iii) H is invariant under the action of the group Sn of all permutations of 
canonical variables q^, pi, i — 1, . . . , N. 

(iv) H admits at least one integral of motion quadratic in momenta and 
functionally independent of H and P. 

To this end we make the following observations: 

a) the hamiltonian (JTJ is time - reversal invariant: qi — > qi, pi — > —pi, i = 
1, . . . , N, is the symmetry of H. Therefore, one can assume from the very 
begining that the integrals of motion are either even or odd functions of 
momenta. In particular, our quadratic integral can be written as 

I(q,p) = I(q,p)+W(q), (3) 
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where I(q, p) is quadratic homogeneous function of momenta. 

b) I(q, p) is itself an integral of motion in the free case (V = 0). This follows 
by considering the highest - order term (in momenta) of /. Now, for a free 
motion any integral is a function of momenta and the combinations (angular 
momentum components) qiPj — qjPi, i, j = 1, . . . , N. Consequently, I is 
also at most quadratic in g's. 

c) Sn~ symmetric quadratic function I(q, p) can be written as a function of 
the variables Ylf=i QdPi \ 0£, = 0, 1, 2. 

Using (a) - (c) one easily finds that the integral / can be written in the 
following general form: 

/ N N N \ 

% p) = « (E a& 2 ) - (E QiPi) 2 + 

V i=i i=i i=i I 

/ N N N \ 

+p(2QP(j: qi P,)-Q 2 (Ep")-P 2 (E^)+ (4) 

V i=l i=l i=l } 

+7 nE^)-^(E^) )+W(q) 

\ i=l i=l I 

where Q = J2?=i Qi- 

Imposing the condition 1 = one arrives at the set of equations determining 
W. The integrability conditions for this set read: 



(D + 2)V = 0; (5) 



here 



9 d 

M kj = q k - qj — 

dqj dq k 

N a 

are the generators of rotations and dilatations, respectively ( j, k = 1, . . . , N 
) 

It is not difficult to find the general solution to eqs. Keeping in mind 

that V is defined up to an additive constant one finds: 
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1) if a 7^ (3N or 7 7^ V is a translationally invariant homogeneous function 
of degree —2. 

2) if a = (3N and 7 = then 



V(q) = V(q) + U £fe-?;) 2 h (7) 




where V is a translationally invariant homogeneous function of degree —2 

while U is an arbitrary differentiable function of one variable. 

Having determined V one can find W. The result follows. 

In the case (1) any quadratic integral is a linear combination of the following 

ones: 

h= (j2q kPl ?JP-2QH 

^ = 2 [f^i) h - [f:^) (8) 

1 3 = 2QP ^£ q k p^ - 2Q 2 H - (jr q 2 ^ P 2 

Note that V enters the above expressions only through H. The integrals 
Ji, J 2 , I3 are, however, not functionally independent. Instead, they obey the 
relation 

I 2 + P 2 I 2 + 2Hh = (9) 

In the case (2), when V is given by eq. (J7J) with nontrivial U, there exists 
one independent quadratic integral: 

/ n \ w 

1 4 = 2 N(£ ll) -Q 2 )H- N(£ <lk Pk ) 2 + 

\ k=l ) k=l 

N N 

+2QP(YlkPk)-P 2 (Yll) (10) 

k=l k=l 

where H is obtained from H through the repacement V — > V. Note that, in 
the case U = 0, Z4 can be written as J 4 = iV/2 + ^3- Therefore in this case 
one can take I\ and I4, as independent quadratic integrals. 
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3 Separation of variables 

Typically quadratic integrals of motion are related to the separation of vari- 
ables by point transformation. Let us introduce the Jacobi coordinates 

2/o = (11) 

1 k 

Vk= , =(fcgjb+i-X)g*)» k = 1, . . . , iV — 1 

fk(k + 1 i=i 



and, as a next step the polar coordinates p, #i, . . . , 9 N _ 2 in the space of 
Jacobi vatiables y±, . . . , yjv-i- Then 

AT-1 -i AT 

^e^4e(*- ?i) 2 (12) 

fc=l ZJV i,j=l 

Moreover, the Hamiltonian for the general potential ((7J) can be rewritten in 
the form 



where F comes from the kinetic part while = p~ 2 G. Inserting (jl3j) into the 
Hamiltonian - Jacobi equation we find that H and P are related (as usual) 
to the separation of time and center - of - mass motion while the separation 
of radial variable gives the additional integral 

C= l -F{e_,p e )+G{9) (14) 

Note that C is dilatation invariant. In fact, one readily checks that 

C=±U (15) 
which explains the origin of our quadratic integral. 

If U = we have two independent quadratic integrals. In order to find the 
second one let us write out the solution to the H-J equation: 

S = ^PQ + Sx(p) + S 2 {6) - Et (16) 



where Si 2 obey: 



i m ^) + Gm = ^7 4 ( 17 ) 



2\dp J ' 2Np 2 4 ' 2N 
The equation for Si can be solved explicitly 



P 2 i 

= ± I dpd 2E - 

Eqs. - (HHJ) imply 



Sl = ± j d j2E---—U (18) 



TV TV/) 2 

as _ h 

aP ~ P 2 - 2NH 



(19) 



However, by standard algorithm j3] |^ is a constant. We conclude that Ii 
arises also from separation procedure for H-J equation. 



4 sL{2, R) symmetry 

The result obtained in the case U = can be also understood from the 
point of view of dynamical sL(2, R) symmetry. Due to the fact that V is 
homogeneous function of degree —2 the three functions H, X = J2k=i Qk an d 
Y — Y<k=i QkPk form the sL(2, R) algebra jl] 

{y, H} = 2H 

{Y, X} = -2X (20) 
{H, X} = -2Y 

One can define the action of this algebra in the space of functions on phase 
space as follows j^j: 

(Af)(q,p) = {A(q,p),f(q,p)} (21) 

where A = X, Y, H. Eq. (}2~T]) defines the representation of the algebra ()20|) . 
In particular, the integrals of motion provide the highest - weight vectors. 
If, in addition, such an integral is a homogeneous function of natural degree 
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(with p (q) having weight +1(— 1)) and a polynomial in p, it defines a finite- 
dimensional representation. From this point of view 1% and I2 are highest - 
weight vectors of doublet and singlet, respectively. It is particularly simple 
to trace the origin of I2] it can be written in terms of sL(2, R) generators 
as 

I 2 = 2XH - Y 2 (22) 

and one immediately recognizes that 1% is simply the Casimir function for 
this algebra. 

To trace the origin of I\ let us note the following. If / is the highest - weight 
vector (J = {J, H} = 0, {Y, /} = nl) then the next - to - highest vector 
J = {X, /} obeys 

j = {J. H} = {{X, I}, H} = -2nl; (23) 

therefore J evolves linearly in time. 

Now, if V is a second integral of degree n' it follows from eq. (|23J) that 
n'I'J — nlJ' is also an integral of motion jH]. Applying this reasoning to the 
triple (H, Y, X) and the doublet (P, Q) one obtains I\. 
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